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ABSTRACT 

Aims. To quantify the transient growth of nonaxisymmetric perturbations in unstratified magnetized and stratified non- 
magnetized rotating linear shear flows in the shearing sheet approximation of accretion disc flows. 

Methods. The Rayleigh quotient in modal approaches for the linearized equations (with time-dependent wavenumber) and the 
amplitudes from direct shearing sheet simulations using a finite difference code are compared. 

Results. Both approaches agree in their predicted growth behavior. The magneto-rotational instability for axisymmetric and 
non-axisymmetric perturbations is shown to have the same dependence of the (instantaneous) growth rate on the wavenumber 
along the magnetic field, but in the nonaxisymmetric case the growth is only transient. However, a meaningful dependence 
of the Rayleigh quotient on the radial wavenumber is obtained. While in the magnetized case the total amplification factor 
can be several orders of magnitude, it is only of order ten or less in the nonmagnetic case. Stratification is shown to have 
a stabilizing effect. In the present case of shearing-periodic boundaries the (local) strato- rotational instability seems to be absent. 
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1. Introduction 

The gas in accretion discs is generally though to be tur- 
' bulent (Shakura & Sunyaev 1973). This allows potential 
' energy to be converted into kinetic energy which can then 
be dissipated (e.g., Frank et al. 1992). This in turn can lead 
I to significant amounts of observable radiation. To sustain 
the turbulence, there has to be some instability. This in- 
stability is now generally thought to be the Balbus-Hawley 
or magneto-rotational instability (hereafter referred to as 
MRI; see Balbus & Hawley 1998 for a review). This in- 
stability is linear and local in that it does not rely on the 
presence of boundaries. It exists already in the axisym- 
metric case in the presence of an external vertical field, 
which makes this instability technically easy to study. 

In the nonaxisymmetric case, an azimuthal pattern 
will be sheared out differentially, i.e. patterns in the in- 
ner parts of the disc are advected faster than in the outer 
parts. After an increasing number of orbits, this causes 
more and more rapid variations in the radial direction. 
Loosely speaking, this makes the radial wavenumber time- 
dependent. 

In the presence of radial boundaries or other radial 
non- uniformities, and certainly also in the presence of non- 
linearity, any radial pattern that was initially harmonic 
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will become anharmonic. This produces a spectrum of 
wavenumbers even for monochromatic initial conditions. 
Mathematically, the stability of such a problem can be 
studied by solving a one-dimensional eigenvalue problem 
subject to radial boundary conditions (Ogilvie & Pringle 
1996). Even though boundary conditions are important 
in that case, the results may still be relevant for driving 
turbulence provided the instability is not limited to the 
vicinity of the boundaries. 

In many numerical simulations of accretion disc turbu- 
lence the shearing sheet (or shearing box) approximation 
has been employed (Hawley et al. 1995, Matsumoto & 
Tajima 1995, Brandenburg et al. 1995). This approxima- 
tion represents the other extreme, where boundaries and 
any non-uniformities are strictly removed. Any instabil- 
ity that survives under these conditions is often referred 
to as a "local" instability, even though its onset proper- 
ties may depend on the system size (as is typical of all 
long- wavelength instabilities). 

For nonaxisymmetric solutions, a purely analytic 
treatment of the shearing sheet model is generally impossi- 
ble, because the solutions exhibit a complicated temporal 
behavior that cannot even be approximated by an expo- 
nential time evolution. In fact, all nonaxisymmetric solu- 
tions only exhibit transient growth whose speed of growth 
depends on the instantaneous wavevector. However, for 
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the MRI it turns out that the Rayleigh quotient ob- 
tained from the time-dependent nonaxisymmetric solution 
is a good approximation to the usual eigenvalue in the 
much simpler axisymmetric problem. The purpose of the 
present paper is to attempt a more systematic survey of 
the nonaxisymmetric problem by studying the dependence 
of Rayleigh quotient in different two-dimensional parame- 
ter planes. We begin with the fairly well understood MRI 
and then turn to the less well understood problem with 
vertical density stratification, but no magnetic field. The 
latter case is less well understood, although it has been 
shown that in the presence of radial boundaries there is a 
linear instability of the Taylor-Couette problem with den- 
sity stratification along the axis (Molemaker et al. 2001, 
Shalybkov & Riidiger 2005, Umurhan 2005). This "strato- 
rotational" instability (hereafter refereed to as SRI) has 
been confirmed numerically in the presence of boundaries 
(Brandenburg & Riidiger 2006), but it may still be local 
in character, i.e. its properties may not be sensitive to the 
presence of boundaries (DubruUe et al. 2005). 

It turns out that in the nonaxisymmetric case without 
boundaries there can at most only be transient growth. 
This is also true of the MRI, which has an instantaneous 
growth rate quite analogous to that in the axisymmetric 
case. In shearing sheet simulations, sustained instability 
can only be the result of nonlinearity allowing mode cou- 
pling and hence the recycling of energy into new modes 
viable of repeated growth. 

2. The shearing sheet formalism 

The full hydrodynamic and magnetohydrodynamic equa- 
tions can always be written in the form 



(1) 



Dq 

Bt 

where q is a state vector combining all components of 
velocity, density, entropy, and the magnetic field, and 
D/Di = d/dt + U ■ 'V is the advective derivative. The 
velocity U can be subdivided into an equilibrium solution 
(the mean flow U) and the departures from the mean flow, 
so 



U = U + u. 



(2) 



In the shearing sheet approximation, the mean flow de- 
pends linearly on the cross-stream coordinate, say x, so 
we assume the mean flow to be 



U = (0,S'a;,0)^ 



(3) 



where S denotes the gradient of the shear flow. In a local 
model of a keplerian disc we have S — — where il is 
the local angular velocity. 

Inserting Eq. in Eq. ^ yields 



(4) 



where the second term on the left hand side has an explicit 
X dependence. It turns out that none of the other terms 



have an explicit x dependence. This is because we are 
restricting ourselves only to linear shear flows, so the flow 
has constant gradients. 

Equation I^J can be linearized with respect to the de- 
partures from the equilibrium solution, q, so we write 
q = q + q' and have 



dq' , r, dq' , 
+ Sx— = Lq , 
at ay 



(5) 



where L is a matrix with differential operators and con- 
stant coefficients. Equation (0) can be solved by mak- 
ing the "shearing sheet" ansatz (Goldreich & Lynden-Bell 
1965, Balbus & Hawley 1992a) 



q'(x, y, z, t) = q[t) exp [ik^{t)x + ikyy + ik^z] . 



(6) 



Note that by differentiating Eq. (|HJ with respect to t, one 
pulls down a term proportional to i{dkx/dt)x. This ex- 
plicitly a;-dependent term can be arranged to cancel the 
second term of Eq. © by choosing 

kx{t) = kxQ - kySt. (7) 
This leads to a set of ordinary differential equations, 

(8) 



dq 

dt ='-'^' 



where L is a matrix with coefficients that are independent 
of X and depend at most only on t. 

In practice, we solve Eq. ((SJ numerically and moni- 
tor the evolution of the norm, {q\q), and of the Rayleigh 
quotient 



m = 



(g|Lg) 



(9) 



N 



where (a|6) = a*bi defines a scalar product, and N is 

the rank of the matrix L. We recall that, if the matrix was 
independent of t, the Rayleigh quotient would be between 
the largest and the smallest eigenvalue of L. 

In the following we discuss first the MRI and turn then 
to the case with stratification and address the possibility 
of the strato-rotational instability (SRI). 



3. MRI 

3.1. Basic equations 

Stratification is unimportant for the MRI, so we focus 
on the simple case with uniform background density. To 
simplify the problem further, we assume an isothermal 
equation of state, so the pressure is given hy p = c^p. 
Here, Cg is the sound speed which is assumed constant. 
The full set of equations, in the presence of rotation with 
angular velocity fl, is then 



pVt/j + J X B, 



(10) 
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~Dt 



B VU BV U, 



Bt ^ 



(11) 
(12) 



where D/Di = d / dt + U ■ 'V is the advective derivative 
with respect to the total flow velocity, B is the magnetic 
field, J = V X i?//io is the current density, and fiQ is the 
vacuum permeability. 

In an accretion disc, tp — ^n^{3x^ — z^) is the tidal po- 
tential that is derived by linearizing the gravitational po- 
tential with respect to some point in the midplane of the 
disc some distance away from the central object. However, 
in the following we ignore vertical gravity and assume a 
more general body force giving rise to the shear flow, so 
we assume ip — flSx^. The equilibrium solution is then 
given hy U = {0,Sx,0), B = const, and p ~ const. The 
mean flow is obtained by balancing 2n x U against V?/'. 
Inserting Eq. ^ into the U ■ VC/ nonlinearity of the mo- 
mentum equation, we obtain 



U ■ VU = U ■ Vu + u ■ VU + u ■ Vu, 



(13) 



where U ■ 'VU — has been used. Note also that the term 
u ■ VJ7 = {0,Sux,0) can be combined with the Coriolis 
force — 2r2 x tt to give the force 



/(«) 



2flUy 

-{2n + s)ux 





(14) 



which describes epicyclic deviations from purely circular 
motion. In terms of the departures from the mean flow, 
u, the momentum equation (|1(J|) can then be written as 
(Brandenburg et al. 1995) 



Dm 



pf{u)-clVp + J X B. 



(15) 



Likewise, in the induction equation, there is the stretch- 
ing term B ■ W on the right hand side of the induction 
equation ifTT)) . which leads to a term B W = (0, SBx,Q). 
We emphasize again that in the two expressions (0, Sux, 0) 
and (0, Sbx, 0), shear only introduces terms with constant 
coefficients. 

So far, all equations have been fully nonlinear. We can 
now linearize the equations about H = 0, B = const and 
denote the departures from the equilibrium solution by a 
prime. The linearized Lorentz force can be written in the 
form 



(V X B') X B = B • VB' - V(B • B') 



(16) 



In the following two subsections we consider the cases of 
imposed fields that point either in the vertical or in the 
azimuthal directions. The former case is particularly in- 
structive, because it allows an instability already in the 
much simpler axisymmetric case. 



3.2. Vertical field and axisymmetric perturbations 

For a vertical field, B = (0, 0, B^), we have 

_ _ 'd,B',-dxB', 

V{B ■ B') = B, I d.B' - dyBi I , (17) 



B VB' 







and the terms on the right hand side of the linearized 
induction equation are 



B ■ Vu' BV u' = B^ 



dzu' 



(18) 



-dxu' 
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With these preparations we can write down the matrix 

L for the MRI with an imposed vertical equilibrium 
field, 

/ 



V 

where we have used the abbreviations 257^ = 2f7 + 5, 
kf- = kiVA (sometimes also ^a), kf = kiCs (or fee), and 
the state vector is q — {vx,Vy,Vz,bx,by,bz,A)'^, where 
the hats denote the shearing sheet expansion analogously 
to Eq. ©. Here we have used the rescaled variables A = 
Cs In p and b — B / ^/JIoPo to make L hermitian if S" = 0. To 
ensure solenoidality of the magnetic field, we calculate bx 
for the initial perturbation from (Balbus & Hawley 1992a) 

bx ^ ~{kyby + kzbz)/kxQ {forkxo^O). (19) 

For kx = ky = 0, the matrix L has an unstable eigen- 
value X = at kz — In Fig. ^ we present a solu- 

tion of the ordinary differential equation (jSj) using L for 
S = — ff^ and VAkz — 1. Note that Re{X/fl) approaches 
|f2, i.e. the most unstable eigenvalue. This suggests that 
the Rayleigh quotient X(t) is indeed a good approxima- 
tion to the most unstable eigenvalue. In agreement with 
earlier work (Balbus & Hawley 1992a, Kim & Ostriker 
2000), the maximum growth rate agrees with the Oort A- 
value (Balbus & Hawley 1992b), which is -^S, or |il for 
keplerian rotation. 



is inde- 



In the special case ky = (axisymmetry), L 
pendent of t and the solution is given by 

q = qe^^ (for fc^, = 0). (20) 

This leads to the well-known dispersion relation (Balbus 
& Hawley 1991) 

[A'^ + X^{2kl + K^) + kiiki ~ 20.3)] (A^ + kl) = 0, (21) 

where = AOQ,^ is the square of the epicyclic frequency. 
The resulting growth rates are 



A| = ~vlk^ - ^ /4^2 fc2j^2 + 1^4 



(22) 



In the range < v\k^ < 2ns, X can have real values, 
where A+ > 0, corresponding to instability. 
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Fig. 1. Axisymnietric MRI: evolution of the real part of 
A (upper panel) and the norm (or gain; see lower 

panel), for v^kz — 1 using a vertical initial field pointing 
in the z direction, with k^ — ky = and S = — |^^- 
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Fig. 2. Nonaxisymmetric MRI: evolution of the real part 
of A (upper panel) and the norm (lower panel), 

for VAky/fl = 1 using an azimuthal initial field pointing 
in the y direction, and kz — 0. 



3.3. Azimuthal field, nonaxisymmetric perturbations 

Next we turn to the nonaxisymmetric problem. Of par- 
ticular interest is the case of a purely azimuthal field 
B = {0, By, 0), which gives rise to the matrix 
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Since L is time dependent, the solution q will not have 

an exponential time dependence and the eigenvalues of 

"(y) 

L cannot be interpreted as growth rate. Therefore we 
use the Rayleigh quotient A instead. 

In Fig. 121 we show the evolution of ReA as a function of 
time, except that time is translated into a corresponding 
variation of kx{t); see Eq. Q. In the second panel we 
show the corresponding variation of (q|q)^/^, where we see 
an increase over about 8 orders of magnitude during the 
time interval in which ReA is systematically positive. Note 
also that max(ReA) « jfl (upper panel of Fig.|2J), which 
is indeed the maximum growth rate of the axisymmetric 
MRI. This supports our interpretation that the Rayleigh 
quotient provides a convenient and quantitative means of 
estimating the growth rate of the instability. 

In Fig. 131 we show a two-dimensional parameter survey 
in the {k,S) plane of the gain factor and of max(ReA). 
Note that both the gain factor and max(ReA) increase 
toward more negative values of S. Remarkable is the fact 
that on the Rayleigh line, S/fl = —2, both quantities vary 
smoothly and do not show any special behavior. 



2.0 




-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 

S / n 



Fig. 3. Nonaxisymmetric MRI: dependence of the gain 
factor in logarithmic representation (upper panel) and of 
maxt(ReA)/r2 (lower panel) on S/fl and VAk/fl for a ver- 
tical field. Note that maxt(ReA)/il reaches the value of 
^S/il for different values of S/Q. 
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Fig. 4. Transient amplification of tlie magnetic field by the 
nonaxisymmetric MRI with an azimuthal field. The solid 
line shows the result from the Rayleigh quotient method 
while the broken lines give the result from direct three- 
dimensional simulations with zero viscosity and zero resis- 
tivity. The square root of the Rayleigh quotient has been 
scaled by a factor Cq = 3.5 x 10^'' to make it overlap with 
i?rms curve. A resolution of only 32'^ meshpoints is com- 
pletely insufficient to resolve even the beginning of the in- 
stability. At least 256^ meshpoints are required to resolve 
the maximum (thick dashed line). After ili > 17 even the 
simulation with 256"^ meshpoints becomes under-resolved. 



3.4. Direct numerical verification 

It is instructive to compare the present modal analysis 
with a direct three-dimensional simulation of the fully 
nonlinear equations in real space. This is done in Fig. ^ 
where we compare the increase of with the re- 

sulting evolution of the root-mean-square magnetic field 
from a direct simulation of the shearing sheet equations. 
Here we also adopt an isothermal gas with constant sound 
speed. 

The initial condition for the 3-dimensional direct sim- 
ulation is obtained by evolving linearized shearing sheet 
equations for ky — 1 and = 10 to the point where 
kx{to) = ^5. [The size of the domain is (27r)'^.] For 
definitiveness, we reproduce here the numerical values in 
Eq. ®: 



-0.311 - 0.037i\ 
-0.461 - 0.054i 
-0.088 -O.OlOi/ 



-0.068 
-0.039 
-1-0.035 



0.581i' 

0.334i ) , (23) 
0.295i, 



and the logarithmic density is given by A = 0.042 — 
0.3647i. The amplitude is chosen to be A = 10~*. FigureEl 
shows images of on the periphery of the simulation do- 
main at different times^. The simulations have been car- 
ried out using the Pencil Code^ which is a high-order 
finite-difference code (sixth order in space and third or- 



^ Animations of the flow can be found 
jhttp : //www.nordita. dk/~brandenb/movies/disc 
|littp : //www. nordita.dk/ software/p encil- code | 



dcr in time) for solving the compressible hydromagnetic 
equations. 

The way how this transient amplification can lead to 
sustained growth is through mode coupling, which is not 
considered in the present analysis. Relevant mode cou- 
plings could come about either through nonuniformities 
in the cross-stream or x direction and through bound- 
ary conditions, or through nonlinearities. In the shearing 
sheet approximation only the latter is a viable possibility, 
and this is probably the mechanism through which the 
early shearing sheet simulations produced sustained tur- 
bulence (e.g. Hawley et al. 1995, Matsumoto & Tajima 
1995, Brandenburg et al. 1995). 

We conclude this section by emphasizing the one- 
to-one correspondence between numerical solutions of 
the shearing sheet equations and the Rayleigh quotient 
method. The nonaxisymmetric case is similar to the ax- 
isymmetric one in that the maximum growth rate is the 
same. This is achieved for — 0, which, in the nonax- 
isymmetric case, can only be the case at one instance in 
time. 

4. Shear flow with stratification 

We now turn to the case with vertical stratification. In 
order to allow for gravity waves, we need to abandon the 
isothermal equation of state and use instead a perfect gas 
law. In terms of specific entropy, s, we can formulate this 
as 



s = Cyhip — Cpln p + So, 



(24) 



where Cp and Cy are the specific heats at constant pressure 
and constant volume, respectively. The specific entropy is 
only defined up to an additive constant, sq. In the ab- 
sence of heating and cooling, the specific entropy remains 
constant on streamlines, i.e. 



Hf.o. 

Bt 



(25) 



Using Eqs H24(l and H25|l we can rewrite the continuity 
equation 1)12(1 in the form 



Dp ^ 



(26) 



where 7 = Cp/c„ is the ratio of specific heats and V-C = 
has been used. We solve Eqs (I25f) and (|26|l together with 
the momentum equation in the form 

Dm 

P^=pfiu)-Vp + pg, (27) 

where f{u) has been defined in Eq. ((14(1 . These equations 
have an isothermal equilibrium solution (denoted by an 
overbar). 



u — 0, s{z) 
where 



H 



(28) 



(29) 
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Fig. 5. Images of the vertical component of the magnetic field, B^, for the nonaxisymmetric MRI with an azimuthal 
field, using different values of t, where the field at i = corresponds to that given by Eq. (|23|l . 



is the pressure scale height and Cg is the sound speed. 
Both H and Cg are constants, and so is po which gives 
the pressure at z = 0. The equilibrium density is given by 
P = IP/cl- 

Linearizing Eqs (|24|I - H27|) about the z dependent equi- 
librium solution yields 



V 



IP 



P ' 



dt 
dp' 



at az 



I dp 



-7pV • It', 



dt 

m-g^^p{^)f{n')-Vp' + p'g, 



(30) 

(31) 
(32) 
(33) 



where primes denote deviations from the equilibrium. 
These four equations have non-constant coefficients. Using 
rescaled variables (denoted by a tilde). 



P = IPope 



-z/2H 



(34) 
(35) 
(36) 



and eliminating p' using Eq. H3()|l . we can rewrite these 
equations in the form 



ds 
'dt 



-r Cs . 



dp 
di 

du 
'dt 



Cs_-. 

dt 7if"^ 



= /(«) 



— + V-U 
2H 



(37) 
(38) 



\2H 



V]p+{p-^^s)^, (39) 



It is now convenient to introduce nondimensional time and 
space coordinates. A natural length scale to choose might 
be H, or better -fH = c^/g. The problem with this is 



that the unstratified limit g — s- is then ill-posed, because 
_ff — > oo . A more flexible alternative is therefore to choose 
a "system" scale L, or better 7L, as our length scale. (The 
length scale L is not to be confused with the matrix L.) 
Thus, we define nondimensional inverse time and length 
scales via 
d Cs d 
dt ~ ^dt 

The final set of equations is then 
dS L 



and V 



(40) 



di 
dp 
'dt 



H 



^7 - 1 



L 

H 



1 - 



du ~ L 
^t^f^'^^-H 



zp — Vp 



1 zs. 



(41) 
(42) 
(43) 



where we have defined f{u) = {2Vl.Uy,—2Q,^Ux,Q)'^ with 
2Vl^ = 2VL + S. Here, h = jLfl/cs and S = -fLS/cs are 
nondimensional angular and shear velocities. We recall 
that g = (0, 0,-17). As a nondimensional measure for the 
degree of stratification we use in the following the symbol 

' ^ - (44) 



9 



H 



The case g = corresponds obviously to the completely 
unstratified case. 

Next, we write these equations in matrix form Q, 
make the ansatz and have a set of ordinary differ- 
ential equations © with the matrix L in the form 

f Q 20. 



-N 
M - \k. 





2f]^ 



\-iK{t) 



20. 









N 





-ifcx(t) \ 

-M - iL 





(45) 



operating on the state vector q — {ux, iiy, u^, s,p)'^ , which 
is defined according to Eq. ©, 



N 



L 
H 



(46) 
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is the nondimensional Brunt- Vaisala frequency, and 

M=^(l-J) (47) 

is another nondimensional number characterizing the de- 
gree of stratification. Note that in the absence of shear, 
5* = 0, the matrix L is skew-hermitian, i.e. Ljj = — L*j, 
and all eigenvalues are purely imaginary or zero. A similar 
nondimensionalization was used by Brandenburg (1988) 
for the case f2 = 0, but space and time coordinates were 
scaled such that the coefficient M became unity. As dis- 
cussed above, this gives an unnecessary restriction in that 
it prevents us from making a continuous transition to the 
unstratified case L/H = 0. Furthermore, in the special 
case 7 = 2, we have M = 0, N = L/H ^ 0, which 
corresponds to the Boussinesq case if we also make the 
assumption of incompressibility. Technically, the latter is 
achieved by multiplying the left hand side of Eq. ^ by 
the matrix diag(l, 1,1,1,0). This means that p can be ex- 
pressed in terms of u, and we are left with only 4 explicitly 
time-dependent equations. 

In the axisymmetric case, ky = 0, or in the absence 
of shear, 5 = 0, the eigenvalues A satisfy the dispersion 
relation XD{X, k) = with 

D^X^ + X^{k^ + K^ + k^) + klN^ + n^{kl + k^), (48) 

where k"^ = k'^ + k'l, and k is the nondimensional epicyclic 
frequency with — AVlVl^ , and k"^ — N"^ -f has been 
introduced for abbreviation. Written in dimensional form, 
this dispersion relation is identical to the usual one for at- 
mospheric waves in an isothermally stratified atmosphere 
(Stein & Leibacher 1974, Ryu & Goodman 1992). There 
are five solutions, one zero-frequency mode, A = 0, and 
four others with ImA = ±a;p(fe) and ImA = ±Wg(fe), cor- 
responding to p- and g-modes, respectively. 

Only in the axisymmetric case, fcy = 0, L becomes 
independent of t [k^ — fc^o = const) and so Eq. H48|l can 
be used even when S ^ Q. From Eq. H48|l it is clear that 
now instability requires 

< ^klN^/ikl + k^), (49) 

which is also known as the Solberg-H0iland criterion (e.g. 
Rudiger & Shalybkov 2002, Narayan et al. 2002). Negative 
values of = 4ilil^ are possible when ~S > 20, corre- 
sponding to Rayleigh's criterion. Instability occurs closest 
to the Rayleigh line when fc^ is small values and k^ large. 

In Fig. El we present the results of a two-dimensional 
parameter survey varying both the degree of stratification, 
g as well as the strength of the shear, S. The Rayleigh sta- 
bility fine, S/fl = —2, is particularly evident in the plot 
of the gain factor (upper panel): for zero stratification the 
gain factor increases sharply as one crosses the Rayleigh 
line. For stronger stratification, the increase in the gain 
factor diminishes, suggesting that stratification has a sta- 
bilizing influence on the Rayleigh- Taylor instability. 

For S/^l > —2, the flow is known to be axisymmetri- 
cally stable. Also in the nonaxisymmetric case the gain 
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Fig. 6. Nonaxisymmetric stratifled case: dependence of 
the gain factor in linear representation (upper panel) and 
maxt(A)/r2 (lower panel) on S/Q and g, for jLkz — 1. 
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Fig. 7. Nonaxisymmetric stratifled case: dependence of 
the gain factor and maxt(A)/r2 on 'jLkxit), for jLkz = 1. 

factor is negligible. Nevertheless, in the (g, S) parame- 
ter plane there are regions with an appreciable maximum 
value of A (lower panel of Fig.|SJ). However, as is seen more 
clearly in the corresponding time trace (Fig.[7|), a positive 
maximum of A is more or less compensated by a corre- 
sponding amount of negative contributions. This explains 
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the absence of a corresponding pattern of max((A) in the 
(g, S) plane. This suggests also that, in contrast to simu- 
lations of the MRI, the transients would be insignificant 
to produce sustained growth. 

5. Conclusions 

The shearing sheet approximation is a powerful tool to 
analyze the behavior of an accretion disc locally. While 
this approximation eliminates a wide class of global insta- 
bihties, it has the advantage of isolating local instabilities 
that are often believed to be responsible for driving tur- 
bulence. Our work has demonstrated that for sufficiently 
weak perturbations there is a direct correspondence be- 
tween the solutions of shearing sheet equations as obtained 
from a numerical code and the Rayleigh quotient obtained 
by solving the modal equations with a time-dependent ra- 
dial wavenumber (t) . 

Within the framework of the shearing sheet approxi- 
mation we have demonstrated the stability of a stratified 
isothermal atmosphere with horizontal keplerian shear to 
fully compressible perturbations. The behavior of the so- 
lution is well characterized by the Rayleigh quotient, X(t), 
which is obtained for a range of different initial conditions. 
The strato- rotational instability (SRI), which has been 
seen in experiments and calculations with radial bound- 
aries (Molemaker et al. 2001, Shalybkov & Riidiger 2005), 
can therefore not be a local instability, as was already dis- 
cussed by Umurhan (2005) and Lesur & Longaretti (2005). 

The present work has also highlighted the absence of 
any correspondence between the growth rates obtained 
from the psewdo-dispersion relation and the actual evolu- 
tion of ReA(t). Of course, there will be transient growth 
from almost any initial condition and for any value of 
S/^l - including S/fl > 0, which would be stable even in 
the presence of magnetic fields (Balbus & Hawley 1991). 
However, under certain conditions this can lead to what is 
known as the bypass transition (Chagelishvili et al. 2003, 
Afshordi et al. 2005), i.e. the transient growth that can 
lead to a new state that is itself unstable. Nevertheless, 
the presence of stratification {g > 0) is here seen to have 
a stabilizing effect, making the connection with the SRI 
implausible. 

We emphasize that in the non-Bous sines q stratified 
case it is important to remove z-dependent coefficients 
before making the exp(ifczz) ansatz. Otherwise, persistent 
real parts of \{t) will arise that refiect merely the fact that 
velocity increases as a wave packet travels into less dense 
regions. 

Finally, it should be pointed out that, if keplerian shear 
flows were locally unstable in stratified media, one might 
have seen this in the fully nonlinear shearing sheet sim- 
ulations of Brandenburg et al. (1995), where the Lorentz 
force was removed and the dynamo-driven turbulence was 
found to decay rapidly; see their Fig. 4. A similar test 
was also done by Hawley et al. (1995), but only for the 
unstratified case which is not relevant here. We note that 
real discs have vertical shear which does allow for a hy- 



drodynamic instability, but it is far less powerful than 
the Balbus-Hawley instability (e.g. Urpin & Brandenburg 
1998, Arlt & Urpin 2004). Other proposals for hydrody- 
namic instability include the baroclinic instability (Klahr 
& Bodenheimer 2003) that has been discussed in the con- 
text of protostellar discs. 
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